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I. INTRODUCTION 

It is often said that the strength of quantum computa- 
tion lies in the phenomena of quantum superposition and 
quantum entanglement. These features of quantum com- 
putation allow to perform the computation on all possible 
inputs that fit to the quantum register. One of the great- 
est achievements in the theory of quantum algorithms is 
quantum search algorithm introduced by Grover. The 
detailed description of this algorithm can be found in [Ij, 

E], E, a- 

Any physical implementation of a quantum computer 
will be error-prone because of the interaction of the com- 
puting device with the environment. In this paper we in- 
vestigate the resilience of Grover 's algorithm in the pres- 
ence of quantum noise. We use the language of density 
matrices and quantum channels. Our goal is to find the 
maximal amount of noise, for which the quantum algo- 
rithm is better, in terms of mean number of operations, 
than classical algorithm. We aim to achieve this objective 
by considering some classes of quantum channels mod- 
elling environmentally induced noise. 

The paper is organised as follows: in Section [IT] we 
make a short review of the subject. In Section [ITT] we de- 
scribe the formalism of quantum information theory. In 
Section IVlwe present the quantum search algorithm. In 
Section Wwe introduce the noise model we have applied 
to the system. In Section |VI| we analyse the results and 
finally in Section [VII| we present some conclusions. 

II. REVIEW OF EXISTING WORK 

The problem of influence of noise on the quantum 
search algorithm has been extensively studied by var- 
ious researchers. In f5| the authors discuss the influ- 
ence of classical field upon a quantum system imple- 
menting Grover's algorithm. The authors of the paper 
[5] ask similar question to the one asked in this work, 
but use Gaussian noise model, which in their case is not 
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described in the language of quantum channels. In |7] 
the authors analyse how imperfections in realizations of 
quantum gates influence the probability of success of the 
quantum search algorithm. In [8J the authors analyse the 
behaviour of the quantum search algorithm realized with 
the use of noisy 7r/4 rotation gates. 

The effect of unitary noise on the quantum search al- 
gorithm is studied in P]. In fTU' the authors examine 
the robustness of Grover's search algorithm to a random 
phase error in the oracle and analyse the complexity of 
the search process. In [IT] the author studies decoherence 
in Grover's quantum search algorithm using a perturba- 
tive method. The authors of [12 use the methods of 
quantum trajectories to study the effects of dissipative 
decoherence on the accuracy of the Grover's quantum 
search algorithm. In [13J the author numerically sim- 
ulates Grover's algorithm introducing random errors of 
two types: one- and two-qubit gate errors and memory 
errors. 



III. FORMALISM OF QUANTUM 
INFORMATION 



A. Dirac notation 

Throughout this paper we use Dirac notation. Symbol 
\'ip) denotes a complex column vector, (^| denotes the 
row vector dual to |?/;). The scalar product of vectors 
I?/;), 10) is denoted by {4>\4>)- The outer product of these 
vectors is denoted as Vectors are labelled in the 

natural way: |0) := (j^), |1) := {\). Notation like \cf)'ip) 
denotes the tensor product of vectors and is equivalent 
to \(h) (g) \ib). 



B. Density operators 

The most general state of a quantum system is de- 
scribed by a density operator. In quantum mechanics 
a density operator p is defined as hermitian {p — ) pos- 
itive semi-definite {p > 0) trace one (tr(p) = 1) operator. 
When a basis is fixed the density operator can be written 
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in the form of a matrix. Diagonal density matrices can 
be identified with probabiHty distributions, therefore this 
formalism is a natural extension of probability theory. 

Density operators are usually called quantum states. 
The set of quantum states is convex [14, and its boundary 
consists of pure states which in matrix terms are rank one 
projectors. Convex combinations of pure states lie inside 
the set and are called mixed states. 



1. Entanglement 

Entanglement is one of the most important phenomena 
in quantum information theory. We say that state p is 
separable iff it can be written in the following form 

M 

1=1 

where > and X]f=i li — ^- ^ state that is not 
separable is called entangled. It is an open problem of 
great importance and under investigation, to decide if 
a given quantum state is entangled or not. 

2. Subsystems 

Given two states p^, of two systems A and B. the 
product state p"^^ of the composed system is obtained 
by taking the Kronecker product of the states i.e. p^^ = 
p"^® p^. 

Let [p"^^\ki be a matrix representing a quantum system 
composed of two subsystems of dimensions M and N . 
We want to index the matrix elements of p using two 
double indices [p'^^lm^ so that Latin indices correspond 

nv 

to the system A and Greek indices correspond to the 
system B. The relation between indices is as follows k = 
(to — 1)7V + /i, I = {n — 1)N + V. The partial trace 
with respect to system B reads tr^fp'^^) — yZ,, P"^f^ = 

^ nfj, 

p^, and the partial trace with respect to system A reads 

Given the state of the composed system p^^ the state 
of subsystems can by found by the means of taking par- 
tial trace of p^^ with respect to one of the subsystems. 
It should noted that tracing-out is not a reversible oper- 
ation, so in a general case 

p^^^tr^(p^^)®trB(p-^^). (2) 

C. Completely positive trace-preserving maps 
(CPTP) 

We say that an operation is physical if it transforms 
density operators into density operators. Additionally we 
assume that physical operations are linear. Therefore an 
operation $(•) to be physical has to fulfil the following 
set of conditions: 



1. For any operator p its image under operation $ 
has to have its trace and positivity preserved i.e. if 
tr(p) = l,p > 0,p = p^ then tr(<I>(p)) = l,$(p) > 
0,$(p) = <f(p)t. 

2. Operator $ has to be linear: 

3. The extension of the operator $ to any larger di- 
mension that acts trivially on the extended system 
has to preserve positivity. This feature is called 
complete positivity. It means that for all positive 
semi-definite p, ^ > the following holds 

($ ® Idim (c ) ) (p ® = * (/o) ® e > 0. (4) 
CPTP maps are often called quantum channels. 

1. Kraus form 

Any operator <J> that is completely positive and trace 
preserving can be expressed in so called Kraus form 
|14] . which consists of the finite set {Ek} of Kraus op- 
erators - matrices that fulfil the completeness relation: 
J2k^k^^k = I- The image of state p under the map <1> 
is given by 

$(p) = ^i?fcpi?fet. (5) 

k 

D. Measurement 

Quantum states cannot be observed directly. In the lit- 
erature one considers two main types of measurements: 
Von Neumann measurement and POVM (Positive Op- 
erator Valued Measure) measurement. In this paper we 
use only Von Neumann measurement but for the sake of 
completeness we also define POVM measurement. 

The mathematical formulation of Von Neumann mea- 
surement is given by a map from a set of projection op- 
erators to real numbers. 

Let us consider an orthogonal complete set of projec- 
tion operators P — {Pi}fLi and the set of real measure- 
ment outcomes O — {oi}fLi. Mapping P O is called 
Von Neumann measurement. Assuming the system is in 
the state p, the probability pi of measuring outcome Oi 
is given by the relation pi = tr(Pip). 

POVM measurement can be considered as a generali- 
sation of Von Neumann measurement. Let us take a set 
of positive operators F — {Fi}fL^ such that X^iLi — ^ 
and the set of real measurement outcomes O — {oi}^i. 
Mapping P — > O is called POVM measurement. Given 
the system is in the state p, the probability pi of measur- 
ing outcome Oi is given by the relation pi — iv{Fip). 
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IV. OVERVIEW OF THE GROVER'S 
ALGORITHM 

Grover's unordered database search algorithm is one of 
the most important quantum algorithms. This is due to 
the fact that many algorithmic problems can be reduced 
to exhaustive search. 

The main idea of the algorithm is to amplify the prob- 
ability of the state which represents the sought element. 
The algorithm is probabilistic and may fail to return the 
proper result. Fortunately the probability of success is 
reasonably high. 



2. Diffusion 

The operator D rotates any state around the state 

i^) = 4^e'i-)' (11) 

V ^ x=Q 

where D can be written as 

D = -i?®"(2|0)(0| = 2|V')(?M - I. (12) 

3. Initialisation 



A. The problem 



Let X be a set and let f : X ^ {0, 1}, such that 



(6) 



for some marked xq € X. 

For the sake of simplicity we assume that X is a set 
of binary strings of length n. Therefore \X\ — 2" and 
/ : {0, 1}" — {0, 1}. We can map the set X to a set 
of states over C*^^ in the natural way: x o \x), form- 
ing orthogonal, complete set of vectors. The goal of the 
algorithm is to find the marked element. 



B. The algorithm 

The Grover's algorithm is composed of two main pro- 
cedures: the oracle and diffusion. 



1 . Oracle 

By an oracle we call a function that marks one defined 
element. In the case of Grover's algorithm, the marking 
of the element is done by the negation of the amplitude 
of the sought state. 

With the use of elementary quantum gates the oracle 
can be constructed using ancilla \q) in the following way: 



0\x)\q) = \x)\q(Sfix)), 



(7) 



where denotes addition modulo 2. If the register \q) is 
prepared in the state 



\q)=H\l) 



|0) 



V2 



(8) 



where H denotes the Hadamard gate, then, by substitu- 
tion, Eq. I?! can be written as 

o\x)\^^i-i)n^W^. (9) 

By tracing out the ancilla we get 

0|x) =-(-l)-^(")|a;). 



(10) 



We begin in the ground state |0 . . . 00) . In the first step 
of the algorithm we apply the Hadamard gate i^"^" on 
the entire register. This transforms the initial state into 
flat superposition of computational base states: 

iJ«^"|0 . . . 0) = ^ (|0 . . . 00) + . . . + |1 . . . 11)) . (13) 



4- Grover iteration 

The core of the algorithm consists of the applications of 
so called Grover iteration gate G = D O. This procedure 
causes the sought state to be amplified and others states 
to be attenuated. 



5. Number of iterations 

The application of the diffusion operator on the base 
state I a;) gives 



D\. 



\xa) 



N 



(14) 



The application of this operator on any state gives 
i^|x)=^a,(-|2:)-^|2/^|y)) 



where 



^(-a, + 2s)|x), 



s = — a,- 



(15) 



fc-fold application of Grover's iteration G on initial 
state \s) leads to |il[T5] 



(16) 



with real coefficients 
1 



COS {2k + 1) e, pk = sin {2k + 1) 0, (17) 
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where 6 is an angle that fulfils the relation 



sin(6l) = 



1 



N 



(18) 



Therefore the coefficients ak , Pk are periodic functions of 
k. After the series of iterations Pk rises. The influence of 
the marked state \xq) on the state of the register results 
in the evolution of the initial state \s) towards the marked 
state. 

The Pk attains its maximum after approximately ^^/N 
steps. The number of steps needed to transfer the initial 
state towards the marked state is of order 0{\/N). In 
the classical case the number of steps is of order 0{N). 



6. Measurement 

The last step of the Grover's algorithm is Von Neu- 
mann measurement. The probability of obtaining the 
proper result is |/3j;p. 



V. NOISE MODEL 

The above discussion of quantum search algorithm has 
been conducted using state vector formalism. In order 
to incorporate the noise into the quantum computation 
model we need to make use of density operators which 
define the quantum state in the most general way. 





' Q -i' 




'1 


(Jy = 


i 


,0-2 = 


-1 



are Pauli matrices. 

Amplitude damping The amplitude damping 
channel transforms |1) into |0) with a given probabil- 
ity a. State |0) remains unchanged. The set of Kraus 
operators is following 



1 







^/r 



a 









Phase damping Phase damping is purely quantum 
phenomenon which describes the loss of quantum infor- 
mation without the loss of energy. It is described by the 
following set of Kraus operators 



1 







a/S 



Bit flip The bit flip family of channels is the quan- 
tum version of classical Binary Symmetric Channel. The 
action of the channel might be interpreted in the follow- 
ing way: it flips the state of a qubit from |0) to |1) and 
from |1) to |0) with probability a. Kraus operators for 
this family of channels consist of a matrix proportional 
to the identity and a matrix proportional to the negation 
gate 



{Vl - al, VofCT^} 



A. Quantum noise 

Microscopic systems that are governed by the laws of 
quantum mechanics are hard to control and at the same 
time, to separate from the environment. The interaction 
with the environment introduces noise into the quantum 
system. Therefore any future quantum computer will 
also be prone to noise. 

a. One-qubit noise There are several one-parameter 
families of one-qubit noisy channels that are typically 
discussed in the literature ^16j. We present them briefly 
below. 

Depolarising channel This is a bi-stochastic chan- 
nel that transforms any state into maximally mixed state 
with a given probability a. The family of channels can 
be deflned using a four-element set of Kraus operators 



where 



" 1 


0' 




■ 


1 " 





1 




1 






Phase flip The phase flip channel acts similarly to 
bit flip channel with the distinction that ctz gate is ap- 
plied randomly to the qubit 

{ \/l — al, y/auz } • 

Bit-phase flip The bit-phase flip channel may be 
considered as joint application of bit and phase flip gates 
on a qubit. Its Kraus operators form is as follows 

{ Vl - Oil, ^foLOy } . 

In all the above families of channels the real parame- 
ter a G [0,1] can be interpreted as the amount of noise 
introduced by the channel. 

h. Multiqubit local channels Our goal is to extend 
the noise acting on distinct qubits to the entire regis- 
ters. We assume that the appearance of an error on a 
given qubit is independent from an error appearing on 
any other qubits. 

In order to apply noise operators to multiple qubits 
we form a new set of Kraus operators acting on a larger 
Hilbert space. 

We assume that we have the set of n one-qubit Kraus 
operators {ek}k^i- We construct the new set of op- 

N 

erators {Ek}k^i that act on Hilbert space of dimension 
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Iterate 



iVf times 



|0) 
|0) 



|0) 



Oracle 



Diffusion 
|0)^|0) 

\x) ^ — \x) 
fora:;>0 



Noise 
Pt + i=*(Pt) 



FIG. 1. The circuit for Grover's algorithm extended with a non-unitary noisy channel. 



2 by applying the following formula 



{Ek} = \J{e^, 



(19) 



where / = {hjl^^ x {i2}l=i x . . . x {iAr}?„=i. 

One should note that the extended channel 
J2k EkpEl is by the definition local 



By applying Eq. (19) to the sets of operators listed 



above we obtain one-parameter families of local noisy 
channels, which we use in further investigations. 



B. Application of noise to the algoritiim 

In order to simulate noisy behaviour of the system im- 
plementing the algorithm we apply a noisy channel after 
every Grover iteration. The evolution of the system is de- 
scribed by the following procedure, which is graphically 
depicted in Fig. [T] 



1. Prepare system in state po 

2. p := i/®"po^f^"^ 



3. Lf vA^J times do: 

(a) apply Grover iteration p := GpG\ 

(b) apply noise p := 

4. Perform orthogonal measurement in computational 
basis. The probability of finding the sought element 
fisp=(^|p|0. 

This approach simplifies the physical reality but it is 
sufficient to study the robustness of the algorithm in the 
presence of noise. In order to study the discussed prob- 
lem we make use of the numerical simulation. Therefore 
some simplification is necessary as the size of the problem 
grows exponentially fast with the number of qubits. 

The tool we use is qusintum-octave |17| . a library that 
contains functions for simulation and analysis of quantum 
processes. 



In our model we assume that it is easy to verify if the 
correctness of the quantum search algorithm the result 
of quantum search results. It is an assumption usually 
made in the complexity analysis of search algorithms. 



VI. ANALYSIS OF THE INFLUENCE OF NOISE 
ON THE EFFICIENCY OF THE ALGORITHM 

An interesting question arises: "What is the maximal 
amount of noise for which Grover's algorithm is more 
efficient than any classical search algorithm?" 

Grover's algorithm is probabilistic, therefore we can- 
not expect to obtain a valid outcome with certainty. We 
assume that if algorithm fails in a given run we will re- 
run it. There is a certain number of reruns for which 
quantum algorithm is worse than classical. We are inter- 
ested only in the statistical behaviour of algorithm and 
calculate the mean value of repetitions. 

Let k = ["y/f be the maximal number of single 
runs of Grover's algorithm for which quantum searching 
is faster than classical. 

We compute Pmin minimal value of success probability 
of single run of Grover's algorithm for which we obtain 
a valid result with confidence C 



Pn 



mm 

p 



{i-{i-pr>c}. 



(20) 



Numerically obtained values ofpmin for confidence level 
C = 0.95 for Grover's algorithm are listed in Tab. |Tj 

For our numerical experiment we assume that sought 
element ^ lies in the "middle" of the space of elements i.e. 
£_ = 2"-i. 

Plots in Fig. [2] depict the influence of noise parameter 
a on a successful run of Grover's algorithm acting on 
six qubits. These values of parameter a for which the 
plots are above threshold level Pmin can be considered as 
the amounts of noise which do not make quantum search 
algorithm less efficient than classical search algorithms. 

We can compare the probabilities from plots in Fig. [2] 
and these for other sizes of quantum registers with Pmin 
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I ' 1 ' 1 ' 1 ' — I I ' 1 ' 1 i ' I 

0.05 0.1 0.15 0.2 0.05 0.1 0.15 0.2 



noise level a noise level a 

FIG. 2. Probabilities of successful run of Grover's algorithm in function of the noise parameter a. Case of six qubits. The 
value for which the plots attain pmin threshold is shown in Tab. |l] 



Size of the system 


k 


Pmin 


TV = 2^ 


1 


0.95000 


N = 2^ 


2 


0.77639 


N = 2"^ 


3 


0.63160 


AT = 2^ 


5 


0.45072 


iV = 2^ 


7 


0.34816 


TV = 2** 


10 


0.25887 



TABLE I. Values of k and Pmin for Grover's algorithm. 



and find the value of the noise parameter a for which 
it is equal to Pniin- The results of the comparison are 
collected in Tab. Ullfor confidence level C = 0.95 and for 
the channels we have described in Section IVl 



C = 0.95 


depolarising 


amplitude damping 


phase damping 


Af = 2* 


0.025 


0.069 


0.177 


Af = 25 


0.032 


0.010 


0.204 


N = 2^ 


0.031 


0.104 


0.190 


Af = 2^ 


0.026 


0.094 


0.158 


Af = 2* 


0.020 


0.075 


0.122 




bit flip 


phase flip 


bit-phase flip 


Af = 2* 


0.025 


0.047 


0.018 


N = 2^ 


0.032 


0.054 


0.024 


N = 2^ 


0.031 


0.050 


0.023 


Af = 2^ 


0.026 


0.041 


0.020 


Af = 2^ 


0.020 


0.031 


0.015 



TABLE IL The maximal values of noise parameter a for which 
Grover's search algorithm is as efficient as classical search 
algorithm in terms of number of uses of the oracle. 

In the case of three qubits we have found that, if we 
expect confidence level C = 0.95 or higher, Grover's algo- 
rithm is never better than classical search algorithm. It 
means that if we want to get the result with high prob- 
ability we need to repeat the quantum search so many 



times that it is more efficient to perform this task classi- 
cally. 

In other cases we have obtained the values of the noise 
parameter a between ~ 0.010 and ~ 0.2 depending on the 
noise type and the size of the system. We observe that 
even if the amount of noise is larger in bigger systems 
(what causes the algorithm to be less efficient) the noise 
is compensated by the quantum speed-up. 

The results gathered in Tab. |TT] do not form a mono- 
tonic pattern. To understand this fact we have to take 
into account that two factors influence these numbers. 
The first one is due to the fact that the same value of 
noise parameter a has larger influence on the quantum 
system for bigger numbers of qubits and for larger N the 
number of Grover iterations and noisy channel applica- 
tions k raises. At the same time the more qubits are used 
to perform the search algorithm the more important the 
quantum speed-up is. 



VII. SUMMARY 

In this work we have shown that a new way of analysing 
the influence of quantum noise on the quantum search al- 
gorithm. Our method uses the model of density matrices 
and quantum channels represented in Kraus form. 

We can conclude that the simulations and analysis have 
shown that only for small amounts of noise the quantum 
search algorithm is still more efficient than any classical 
algorithm. 

From our numerical results we conclude that differ- 
ent forms of noise have different impact on the efficiency 
of the quantum search algorithm. The least destructive 
form of noise is phase damping, more destructive is am- 
plitude damping and the most destructive is the depolar- 
izing channel. 

Further work would have to take into account quan- 
tum error correcting codes and more precise noise models 
dependent on the implementation. One of the research 
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directions would be to analyse the quantum search al- 
gorithm in the framework of control Hamiltonians taking 
into account Markovian approximation of quantum noise. 
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